A large-deviation approach to space-time chaos 
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In this Letter we show that the analysis of Lyapunov-exponents fluctuations contributes to deepen 
our understanding of high-dimensional chaos. This is achieved by introducing a Gaussian approxi- 
mation for the large deviation function that quantifies the fluctuation probability. More precisely, a 
diffusion matrix D (a dynamical invariant itself) is measured and analysed in terms of its principal 
components. The application of this method to three (conservative, as well as dissipative) models, 
allows: (i) quantifying the strength of the effective interactions among the different degrees of free- 
dom; (ii) unveiling microscopic constraints such as those associated to a symplectic structure; (iii) 
checking the hyperbolicity of the dynamics. 

PACS numbers: 05.45.-a, 05.40.-a, 05.10.Gg, 05.45.Jn 



Introduction - There are two complementary reasons 
to investigate the links between statistical mechanics 
and space-time chaos. On the one hand, (equilibrium) 
statistical mechanics provides an effective framework to 
describe the evolution of nonlinear systems. This is 
achieved through the introduction of the so-called ther- 
modynamic formalism [l| and is based on a suitable parti- 
tion of the phase-space and the consequent interpretation 
of the time-axis as an additional spatial direction. This 
approach proved to be very effective in the characteri- 
zation of low-dimensional systems and has contributed 
to establish, e.g., the relationship between Lyapunov ex- 
ponents on the one side and fractal dimension or the 
Kolmogorov-Sinai entropy, on the other [2(]. A general- 
ization of the approach to spatially extended systems is 
formally possible, but almost unfeasible, because of the 
difficulty to construct appropriate phase-space partitions 
Q. On the other hand, a detailed understanding of high- 
dimensional chaos can help bridging the gap between mi- 
croscopic and macroscopic evolution, thereby laying the 
foundations for a dynamical theory of (non) equilibrium 
statistical mechanics. In this perspective, the estima- 
tion of a suitable large deviation function appears to be 
the most promising strategy. This idea proved already 
fruitful in the context of a stochastic dynamics, where 
some exact calculations have been performed in simple 
but non-trivial models of interacting particles HH- In 
the context of chaotic systems, instead, this approach is 
the core of the Gallavotti-Cohen fluctuation theorem Q , 
that is proved under the hypothesis that, in the thermo- 
dynamic limit, the evolution of typical dynamical systems 
is effectively hyperbolic. 

In this Letter, we propose an approach that can con- 
tribute to make progress along both directions, without 
introducing any assumption of the underlying dynam- 
ics. More precisely, we suggest to study the fluctuations 



of the Lyapunov exponents (LEs) along the lines of the 
multifractal theory Q. One of the advantages of deal- 
ing with LEs and their fluctuations (in the long-time 
limit) is that they are dynamical invariants, i.e. they are 
independent of the parametrization of the phase space. 
An exact implementation in generic nonlinear models is 
out of question. Nevertheless, here we show that use- 
ful information can be extracted by working within the 
Gaussian approximation. For instance, we show that the 
(cross)correlations among all pairs of LEs and, in partic- 
ular, their scaling behavior with the system size allows 
estimating the strength of the effective interactions that 
spontaneously emerge among the various degrees of free- 
dom. Notice that our analysis goes beyond the usual ex- 
tensivity assessment of space-time chaos, that is linked to 
the existence of a limit Lyapunov spectrum. In fact, we 
will see that the fluctuations of a chain of contiguous non- 
interacting systems are substantially different from those 
of a typical chain of interacting systems. Finally, our ap- 
proach allows testing the hyperbolicity of the underlying 
dynamics, by: (i) comparing the results obtained for dif- 
ferent dcfintions of the Lyapunov exponents, (ii) testing 
phenomena like the dominance of Oseledec splitting [7[, 
and (iii) quantifying dimension variability Q. 

Theory - Let A, (t) denote the ith expansion factor over 
a time r in tangent space. The rate A; = Aj (r) /r is the 
so-called finite-time Lyapunov exponent (FTLE), which, 
in the infinite-time limit, converges to the LE A^ (here 
and in the following, overlines denote time averages). For 
finite t, FTLEs fluctuate around the asymptotic values. 
The theory of large deviations suggests that, in the long- 
time limit, the probability distribution P(X, r) (where 
A = {Ai , A 2 , . . . , Xn} and N is the number of degrees of 
freedom) scales with r as 



P(X, 



(1) 
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where 5(A) is the positive-definite large deviation func- 
tion whose minimum (equal to zero) is achieved in corre- 
spondence of the LEs Aj . 5 has been mostly studied in 
contexts where A reduces to a scalar variable A, as it hap- 
pens for low-dimensional chaos, where it is known that 
5(A) is itself a dynamical invariant Q. This is because 
FTLE fluctuations originate from passages of a chaotic 
trajectory in the vicinity of periodic orbits with different 
stability properties. There is no reason to doubt that 
dynamical invariance is lost upon increasing the dimen- 
sionality of the phase-space. 

If a system is the Cartesian product of uncoupled vari- 
ables, 5 is the sum of functions, each dependent on a 
single Xi, but interactions bring new terms. Although 
determining 5 is too ambitious a task, relevant features 
can be uncovered by expanding it around the minimum 
A, and retaining the first (quadratic) non-zero term, 



5(A) w -(A-A)Q(A- A)* 



(2) 



(f denotes the transpose), an approximation that is 
equivalent to assuming a Gaussian distribution. In prac- 
tice, it is preferable to consider the symmetric matrix 
D = Q _1 . In fact, the elements Dij can be directly de- 
termined by estimating the (linear) growth rate of the 
(co) variances of (Aj(r) — Ajr), 



= lim 



Ai{r)Aj(T)- \X 3 t' )/t. 



(3) 



There are three basic definitions of FTLEs. One can com- 
pute them:(i) by repeatedly applying the Gram-Schmidt 
orthogonalization procedure to a set of linearly inde- 
pendent perturbations (backward or Gram-Schmidt Lya- 
punov vectors); (ii) by performing this procedure along 
the negative time axis (forward Lyapunov vectors); (Hi) 
by making reference to the covariant Lyapunov vectors 
[9(. In the infinite-time limit, the three methods pro- 
duce identical LEs. For finite but long times, as long 
as FTLE fluctuations are connected to "visits" of differ- 
ent periodic orbits, the three defintions should be again 
equivalent. In fact, we have systematically verified that 
the correlations computed with the methods (i) and (Hi) 
are basically indistinguishable as soon as the diffusive 
asymptotic behavior sets in ll|. This result provides a 



first evidence of the effective hypcrbolicity of the under- 
lying dynamics. In fact, it implies that the differences 
induced by the presence of homoclinic tangencies are so 
rare that they do not affect our perturbative analysis. 

The information contained in D can be expressed in 
a compact form by determining its (positive) eigenvalues 
Hk (k — 1,N), which represent the fluctuation ampli- 
tudes along the most prominent directions. Some of the 
eigenvalues may turn out to vanish because of more or 
less hidden constraints. For instance, in the presence of 
a constant phase-space contraction rate, Xi = const 



and all {A^} n-tuples lie in a same hyperplane. As a re- 
sult, one eigenvalue of D is equal to zero: its correspond- 
ing eigenvector is perpendicular to the hyperplane itself. 
Another instructive case is that of symplectic dynamics: 
since the LEs come in pairs whose sum is zero, the fluc- 
tuations of the negative LEs are perfectly anticorrelated 
with those of the positive ones, so that Dij has an addi- 



tional symmetry S, i.e., D 



D 



i,N+l—j 



= -D t 



Altogether, the possible existence of zero eigenvalues re- 
inforces the choice of studying D rather than its ill- 
defined inverse Q. Moreover, since the matrices Q and 
D are diagonal in the same basis, and the eigvenvalues 
of Q are the inverse of those of D, we can infer the scal- 
ing behavior of the former ones from that of the latter. 
One must simply be careful and discard the redundant 
variables, associated to the zero eigenvalues. In particu- 
lar, since, as we shall see, pk cx 1/N, the large deviation 
function 5 turns out to be proportional to the number 
N of degrees of freedom, i.e., it is an extensive quantity. 

Model analysis - We start the numerical analysis by 
studying a chain of Henon maps [lOj 

x n (t + 1) = a - [x n (t) + eVx n (t)] 2 + bx n (t - 1) , (4) 

where T>x n = (x„_i — 2x n + x n+ \) is the discrete Lapla- 
cian operator. We have chosen a = 1.4, b = 0.3, 
e = 0.025 and used periodic boundary conditions (the 
same conditions have been chosen in the other models 
too). The results are shown in Fig. [TJ In panel a) we 
report the self-diffusion coefficients Da (see the sym- 
bols). The clean overlap of the scaled curves indicates 
that Du(p) ps 1/A 85 . This means that the LEs self- 
average in the thermodynamic limit. The singular be- 
havior exhibited by Dn(p) for p — > 1 follows from the 
different scaling behavior of the first and Ath exponent 
which decrease as l/\/~N. In Fig. [TJd we plot D^ along 
the column j = 2 A/5. The off-diagonal terms decrease 
as 1/N, so that the matrix D becomes increasingly diag- 
onal in the thermodynamic limit. Finally, the eigenvalue 
spectrum pk (see Fig. [lb) decreases like 1/N. This im- 
plies that the eigenvalues of Q are proportional to N, 
i.e. the large deviation function is an extensive observ- 
able. Moreover, the 1/a singularity at a = means that 
the leading eigenvalue p\ does not decrease, i.e., there 
exists one direction in phase-space along which fluctua- 
tions survive even in the thermodynamic limit. The phys- 
ical meaning of this feature is to be understood. Finally, 
the eigenvalue spectrum exhibits a remarkable and unex- 
pected property: half of it is equal to zero. A close in- 
spection of the whole correlation matrix reveals that this 
is because D is 5-symmctric. By further investigating 
the Jacobian matrix J, we have discovered that it indeed 
satisfies the symplectic-like condition JAJ T = — bA (sec 
12 1). Unlike the similar case studied in Ref. [l3j |. here 
A is a generic antisymmetric matrix depending on t. Al- 
together, these results indicate that LEs come in pairs, 
such that Xi + Ajv+i-i = hi b. 
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FIG. 1: (Color online) Diffusion coefficients in a chain of 
Henon (symbols) and symplectic (lines) maps. In all pan- 
els, squares, diamonds, plusses and circles refer to N — 40, 
80, 160, and 320 respectively, while dotted, dashed and 
solid lines correspond to N — 32, 64, and 128, respectively. 
The results have been obtained by iterating the chain over 
5 x 10 6 time steps. Panel a) contains the diagonal elements, 
p — (i — 1/2) /N; panel b) refers to the column j = 2N/5, 
p — i/N; panel c) refers to the eigenvalues of D, ordered from 
the largest to the smallest one, a = k/N. 



Next we have studied a chain of symplectic maps, 

p n (t + l) = p n (t) + K (sin A6 n+1 (t) - sin A9 n (t)) 
0„(t + l) = 6 n {t) + Pn (t + 1) (5) 

where both p n and 8 n are defined modulus 2ir and A9 n = 
n — The model has been simulated for K = 4. 

In Fig. Q] (see the lines), one can notice that the overall 
scenario is very similar to that one observed in the chain 
of Henon maps, including the behavior of the diagonal 
elements. The major difference concerns Dnn which, 
instead of decreasing faster, it now decreases slower than 
in the bulk. 

Finally, we have considered a chain of Stuart-Landau 
oscillators as an example of a continuous-time dissipa- 
tive system. The model can be viewed as the spatial 
discretization of a complex Ginzburg-Landau equation, 
a prototypical model of space-time chaos. The evolution 
equation writes 



(1 + ic)\a n \ 2 a n + (1 + ib)h~ 2 Va n . 



(6) 



We have fixed c = 3, b — —2 and h = 1/2, which corre- 
sponds to a regime of amplitude turbulence [3] ■ In this 
model we cannot draw clear conclusions on the scaling 
behavior of the D elements, because of larger finite-size 
corrections (see Fig. However, the eigenvalues behave 
quite similarly to the two previous cases: (i) the over- 
all spectrum scales as 1/N; (ii) the maximum eigenvalue 
remains finite for increasing N; (iii) a large fraction of 
the spectrum is nearly equal to zero. In this case, the 
singularity is due to the appearancejheyond a certain 
p-value) of pairs of degenerate LEs [14| which fluctuate 



synchronously. Notice also the drops of the diffusion co- 
efficient indicated by arrows 1 and 2 in Fig. [2^, that are 
discussed below. 




FIG. 2: (Color online) Diffusion coefficients in a chain of 
Stuart-Landau oscillators. In all panels circles and squares 
refer to N = 64 and 128, respectively. The panels are or- 
ganized in the same way as in the first figure. The column 
reported in panel b) corresponds to j — N/2. 



Discussion - The common property exhibited by all 
of the three models is the 1/N scaling of the eigenvalue 
spectrum of the diffusion matrix D. This implies that the 
matrix Q and, thereby, the large deviation function S are 
proportional to the number of degrees of freedom, i.e., S 
is an extensive observable. It is interesting to notice that 
such a property holds in spite of the long-range correla- 
tions that are revealed by the strength of the off-diagonal 
terms of D (in all models, they provide a substantial con- 
tribution to the scaling behavior of the eigenvalues) . 

Next we discuss some physical implications of the 
structure of the large deviation function S. We start 
from the occurrence of occasional changes in the order 
of the FTLEs, a feature that is related to the concept of 
dominated Oseledec splitting The splitting is dom- 
inated with index i if there exists a finite tq such that 
Xi > Ai+i for all t > tq [15]. This property implies 
the absence of tangencies between the corresponding Os- 
eledec subspaces identified by the ith vector and the sub- 
sequent one [7] • The probability of order-exchanges can 
be inferred from the fluctuations of SXi = Xi — A^+i. Its 
diffusion coefficient Ki can be expressed in terms of the 
D elements, Ki — Da + D i+ i^ +1 — 2Dj ) ,- + i. Since the 
probability of an exchange of FTLEs is equal to the prob- 
ability P™ 1 of observing a negative SXi, we have, in the 
Gaussian approximation, 



P" 



1 



erfc <^ (Aj - A i+ i] 



T 

2K~ 



(7) 



where "erfc" is the complementary error function. The 
analysis of the three models reveals that, in the bulk, Ki 
scales always as 1/N (see Fig. [3]) [l^|. Since the distance 
between consecutive LEs scales also as 1/N (this follows 



4 



from the very existence of a limit LE spectrum) , one can 
conclude that erfc(ci y / r/N)/2 for some c, > 0. 

This means that in the large N limit, order exchanges oc- 
cur with a finite probability and no dominated splitting 
is present. However, in the Henon maps, at p = 1, there 
is a gap in the Lyapunov spectrum. Therefore, since Kn 
vanishes (as 1/y/W), the probability of order exchanges 
goes to zero, indicating that stable and unstable mani- 
fold are mutually transversal and the system effectively 
hyperbolic. The absence of a gap in the Lyapunov spec- 
trum of the symplectic maps prevents us from drawing a 
similar conclusion in that model. In the Stuart-Landau 
chain, P™ vanishes close to arrow "2" (see Figs. [5^ and 
\Sjp), since = (and A, g£ Aj+i), thus implying that the 
splitting is dominated [l8j- This suggests the existence 
of two transversal subspaces, consistently with the claim 
that the attractor is embedded in a supporting manifold 
containing the physical modes 14 , lij] . Since the dimen- 
sion of the supporting manifold is even larger than the 
Kaplan- Yorke dimension (equal to 0.27Y), we must con- 
clude that the overall dynamics is not hyperbolic. 
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FIG. 3: (Color online) Rescaled diffusion coefficients Ki of the 
LE differences 8X in Henon (symbols) and symplectic (lines) 
maps (panel a) and Stuart-Landau oscillators (panel b), p = 
i/N. Same notations as in the two previous curves. 



Now, we analyse the invariant measure, introducing 
the expansion rate C(T>) = JY Xi of a generic volume 
of dimension T>, over a time r. The Kaplan- Yorke di 
mension T>f is obtained by imposing C(Df) 



[17 1 



Under the assumption of small fluctuations, one can 
express the diffusion coefficient Af of V in terms of 
the analogous coefficient A £ = JY j ^ij °f A by lin- 
earizing the function C(T>) around 2?/. This leads to 
Af = A c /(\ Vf ) 2 - The dimension fluctuations Af, 
can now be estimated by invoking an Ansatz similar to 
Eq. |T]) which, in the Gaussian approximation, writes 
P(D) cx £ (' D - l:) /) 2 /(2Af ) ) w h ere the box-size e must be 
linked to the time variable. By following Ref. (l9j |. it is 
natural to assume that e ~ exp(— |A-p / |r), thereby ob- 
taining Af = Af\r> f = A c /\x> f . In the chain of Henon 
maps, Af rj O.I27V, i.e. dimension fluctuations are ex- 
tensive. This implies that the naive idea is wrong and 
it is necessary to build a more refined picture to refer to 



high-dimensional chaotic attractors. 

Conclusions We have shown that a fluctuation analysis 
can deepen our understanding of high-dimensional chaos. 
The main result is the discovery of a subtle form of exten- 
sivity, i.e. the proportionality of the large deviation func- 
tion to the system size. This result is nontrivial, since it 
arises in a context of effective long-range correlations and 
there are even examples of stochastic models, where the 
large deviation function is not extensive As for the 
discrepancy between the scaling exponent of the diagonal 
elements and of the eigenvalues of D (0.85 vs. I) it is nec- 
essary to study larger sizes to decide whether it is due to 
finite-size corrections. Moreover, our approach provides 
a new way of investigating the hyperbolicity of a given 
dynamics (including dimension variability), although we 
are aware that the last word can be said only by go- 
ing beyond the perturbative approach described in this 
Letter. The method introduced in Ref. 20] to identify 
trajectories with unprobable stability properties, makes 
this perspective not so remote. 
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